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Abstract. We study the combination of hyperfine and Zeeman structure in the spin-orbit 
coupled A'E+ — lr'Y\ u complex of S7 Rbi. For this purpose, absorption spectroscopy at a 
magnetic field around B = 1000 G is carried out. We drive optical dipole transitions from 
the lowest rotational state of an ultracold Feshbach molecule to various vibrational levels with 
0+ symmetry of the A — b complex. In contrast to previous measurements with rotationally 
excited alkali-dimers, we do not observe equal spacings of the hyperfine levels. In addition, 
the spectra vary substantially for different vibrational quantum numbers, and exhibit large 
splittings of up to 160 MHz, unexpected for 0 + states. The level structure is explained to be 
a result of the repulsion between the states 0 + and 0 of b^Yl u , coupled via hyperfine and 
Zeeman interactions. In general, 0 and 0 + have a spin-orbit induced energy spacing A, that 
is different for the individual vibrational states. From each measured spectrum we are able 
to extract A, which otherwise is not easily accessible in conventional spectroscopy schemes. 
We obtain values of A in the range of ±100 GHz which can be described by coupled channel 
calculations if a spin-orbit coupling is introduced that is different for 0 and 0 + of /? , IT„. 
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1. Introduction 

The strongly spin-orbit coupled A 1 L+ — b 3 H u complex of alkali-metal dimers has been studied 
in great detail in recent years, stimulated by the fruitful combination of high-resolution 
spectroscopy and numerical close-coupled calculations. Various homonuclear (Rb 2 [1, 2, 3], 
Cs 2 [4, 5], Na 2 [6, 7], K 2 [8, 9, 10, 11], Li 2 [12, 13]) and heteronuclear (NaRb [14, 15], RbCs 
[16, 17, 18], KRb [19], NaCs [20], KCs [21, 22, 23], NaK [24, 25, 26]) species have been 
investigated and modeled. Potential energy curves as well as r-dependent spin-orbit-coupling 
functions were extracted, where r is the internuclear separation. Concerning the hyperfine 
structure of the A — b state, however, only little experimental data is available so far. 

For thermal and thus rotationally excited samples of Na 2 and K 2 hyperfine structures with 
line splittings up to hundreds of MHz, characterized by nearly equidistant separations of the 
energy levels were observed [7, 11]. Such hyperfine structures of the Q. = 0 components of the 
A — b complex come about owing to the molecular rotation that mixes different Q. components. 
For the case of low rotational angular momentum J, line splittings of at most a few MHz are 
expected. Indications of such small hyperfine splittings for J — 1 RbCs molecules in state 
Q — 0 were reported in Ref. [16], but a detailed analysis was not given. 

In this work, we investigate the combined hyperfine and Zeeman pattern of the A — b 
complex for Rb 2 molecules with J = 1 observed by exciting an appropriate Feshbach 
molecular state [see level scheme in Fig. 1(a)]. Particularly for states, where the main 
component exhibits b 3 n„ 0 + symmetry, we measure large level spacings of up to 160 MHz. 
Furthermore, the line pattern is not equally spaced and the overall structure changes strongly 
from one vibrational level to another. Consequently, our spectra are dominated by a 
mechanism different from the one discussed previously in the context of fast rotating 
molecules. In fact, we find that the observed energy level structures corresponding to 
vibrational states of b 3 TI u 0 + arise from second order hyperfine and Zeeman interaction 
coupling the 0 + and 0 components of b 3 n u . More precisely, these two interactions work 
together in a cooperative way enhancing the effect. By fitting a relatively simple model 
to our data we extract the initially unknown frequency spacing A between 0 and 0 + for 
each vibrational level. This is an important result of our work because the state b 3 U u 0 
is not directly accessible in spectroscopy schemes starting from any singlet or triplet ground 
state molecular level. Our derived values for A systematically deviate by about 90 GHz from 
predictions of close-coupled channel calculations. We interpret this as a difference in the 
spin-orbit coupling function for 0 and 0 + . 

This article is organized as follows. In section II, we give an overview of the experimental 
setup and the spectroscopy scheme. Then, section III describes the relevant molecular energy 
states needed for the presentation of our experimental results in section IV. In section V we 
introduce a simple model that fully explains the characteristics of the observed spectra. Our 
model calculations are discussed in section VI along with the interpretation of the data and 
the determination of A for the investigated vibrational states of b 3 Il u . Finally, in section VII 
we describe the extension of the potential scheme needed for modeling the observations by 
coupled channel calculations. 
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Figure 1. (a) Spectroscopy scheme. Weakly bound Feshbach molecules are irradiated by a 
laser pulse and excited to molecular levels of the A 1 — b 3 Yl u manifold from where they 
spontaneously decay to nonobserved states. The potential (2) 3 £+ is included because it 
couples to 0 of b 3 n u (see text). Furthermore, the inset shows the level structure in the 
vicinity of the Feshbach resonance (FR). At a magnetic field of B = 999.9 G the Feshbach state 
(indicated by the black circle) is located 1.748GHz x h below the | f a = 1 ,m.f a = 1) + \fb = 
1 ,rrif b = 1) dissociation threshold at 0G. In (b), the vibrational levels v A and v/, within the 
A*£+ — b 3 n„ 0+ complex that are relevant for our measurements are depicted. All potential 
curves are taken from [27], while the energies of v^ and v/, correspond to the calculated values 
given in [28] (see also tables 1 and 2). 
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2. Experimental setup 

We carry out our experiments with a pure sample of about Nq = 3 x 10 4 weakly bound 87 Rb 2 
Feshbach molecules which have both A 'Ej and a 3 L + character. The setup and the molecule 
preparation scheme are described in detail in Refs. [29, 30]. Therefore, they are just briefly 
presented here. Initially a BEC or ultracold thermal cloud of spin-polarized 87 Rb atoms with 
total angular momentum / = 1, mj — 1 is loaded into a rectangular, 3D optical lattice at a 
wavelength of A = 1064.5 nm, which is formed by a superposition of three linearly polarized 
standing light waves with polarizations orthogonal to each other. By slowly crossing the 
magnetic Feshbach resonance (FR) at 1007.4 G from high to low fields, pairs of atoms in 
doubly occupied lattice sites are converted into weakly bound molecules. Afterwards, the 
magnetic field is set to 999.9 ± 0.1 G, where we perform the spectroscopy. In order to get rid 
of remaining atoms, a combined microwave and light pulse is applied which removes them 
from the lattice. We end up with a pure ensemble of molecules that resides in the lowest 
Bloch band of the optical lattice with no more than a single dimer per lattice site. The lattice 
depth for the molecules with respect to each of the standing light waves of the optical lattice 
is about 64 Er where Er = hr/ (2mA 2 ) represents the recoil energy. Here, m denotes the mass 
of the molecule and h is Planck’s constant. Since at these lattice depths the tunneling rate is 
very small, intermolecular collisions are strongly suppressed, and we measure lifetimes on 
the order of 1 s. 

Figure 1(a) shows the spectroscopy scheme. The Feshbach molecule ensemble is 
irradiated by a rectangular light pulse for a duration T of typically a few ms. At the location of 
the molecular sample the beam waist is about 1.1 mm. For the observed spectra, we used laser 
powers of tens or hundreds of gW. The light propagates orthogonally to the quantization axis 
which is defined by the applied magnetic field that points in vertical direction. By using a 
half-wave plate we can choose the light being polarized either in the horizontal plane or in the 
vertical axis giving rise to a transitions (i.e., (J 1 and c> ) or K transitions. Molecules, that are 
resonantly excited from the Feshbach state to a level of the A—b complex, are in general lost 
due to subsequent fast decay to nonobserved states. We measure the remaining fraction N/Nq 
of Feshbach dimers. For this purpose, we dissociate the molecules by ramping back over the 
Feshbach resonance and detect the corresponding atom number via absorption imaging. 

The spectroscopy is performed at wavelengths between 1042 and 1068 nm (correspond¬ 
ing to about 9360 — 9600 cm -1 ) using a grating-stabilized cw diode laser that has a short-term 
linewidth of ~ 100kHz. This laser is frequency-stabilized to a Fizeau interferometer waveme- 
ter (High Finesse WS7), with an update rate of about 10 Hz. As the laser frequency drifts 
between updates, we obtain a frequency stability of ±(2 — 5)MHz. The wavemeter is cali¬ 
brated to an atomic 87 Rb reference signal at 780 nm in intervals of minutes. It has a specified 
absolute accuracy of 60 MHz, but the accuracy is on the MHz level for difference frequency 
determinations within several hundred MHz. Furthermore, over a period of several months 
we checked the frequency readings of the wavemeter for the same molecular transitions and 
did not find deviations of more than ±10 MHz. This demonstrates the good reproducibility of 
the wavemeter readings in connection with the calibration mentioned above. 
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3. Relevant states 

3.1. Feshbach molecules 

The Rb 2 Feshbach molecules in our experiment are weakly bound dimers with both singlet 
and triplet character, i.e., the selected state is a mixture of A ] E+ and <r/ 3 E+ ([30, 31]). 
However, only the A ! E+ component allows to drive transitions to the A —b complex because 
for an electric dipole transition the u/g symmetry has to change and the A — b complex has 
u symmetry. According to coupled channel calculations, at a magnetic field of 999.9 G the 
singlet component, mainly characterized by S = L = R = 0, 1 = 2, mi = 2, F — 2, contributes 
16% to the Feshbach state which has the exact quantum numbers mp = 2 and parity +. Here, 
S, L, R. 1 and F (F = R+L + S +1) denote the quantum numbers of the total electronic spin, 
the total orbital angular momentum, the rotation of the atom pair, the total nuclear spin, and 
the total molecular angular momentum, respectively. Furthermore, m/ and nip represent the 
corresponding projections onto the quantization axis. Consequently, the singlet component of 
the Feshbach molecules has J = 0 (J = R+L + S). 

The inset of Fig. 1(a) shows the molecular level structure in the vicinity of the Feshbach 
resonance. Throughout the present work, all excitation energies are given with respect to the 
| f a = 1 ,nif a = 1) + \ fh — 1 ,nif h = 1) atomic dissociation limit at 0G. Note, its energy is 
8.543 GHz x h below the atomic dissociation limit when hyperfine interaction is ignored. At 
a magnetic field of 999.9 G the Feshbach state is located at —1.748 GHz x h. Here, the main 
contribution is determined by the Zeeman shift of the atom pair | f a = 1, m/ a — 1) + //, = 
1, mj h = 1). The molecular binding energy is only about 20MHz x h with respect to this 
threshold. 

3.2. A 1 E+ — b 3 II u complex 

Spin-orbit interaction leads to a mixing of the states A ! E+ and b 3 U u forming the A — b 
complex. In a simple approach this mixing comes about in two steps. First, due to spin-orbit 
coupling the state b 3 H u splits up into three components, Q = 0,1,2. The quantum number 
Q. denotes the projection of the sum of all electronic angular momenta onto the intemuclear 
axis and equals the projection of the molecular angular momentum J on the same axis. For 
Rb 2 the relative separation of the three terms is about 80 cm as mainly determined by 
the atomic spin-orbit splitting of Rb in its 5 2 P state. At this stage, the b 3 Yl u , £1 = 0 state 
has two degenerate components, 0 + and 0 _ . Second, spin-orbit coupling mixes A 'E+ (i.e., 
0 + symmetry) and b 3 H u 0 + , whereas the b 3 Tl u 0~ component couples to (2) 3 E+ 0~ [see 
Fig. 1(a)]. As a consequence of the repulsive interactions 0 + and 0“ of /? 3 n„ are separated 
from each other, which is referred to as A-type splitting [32]. This effect is crucial for the 
interpretation of the observations of the present work. 

The vibrational levels of the A — b states relevant to our measurements are illustrated 
in Fig. 1(b). The levels with dominant triplet (singlet) character are indicated by vibrational 
quantum numbers v/, (v^). Moreover, tables 1 and 2 list the numerical values for the term 
energies and the b state admixtures calculated by Drozdova el al. [1] and taken from [28]. The 
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Table 1. Comparison of calculated (Scale) and measured (S e xp) level energies for various 
vibrational levels v ;1 of the A'E+ state with / = 1. All level energies S exp are observed with 
^■-polarized light. The column e = (S ca i c — E exp )/hc gives the difference of the measured and 
predicted values. Furthermore, the parameter pi, denotes the admixture of the Ir' 1 FT,, potential 
and 8 represents the measured frequency difference between the <y and the 7Z resonance. For 
the case of = 67, 68 and 70 we only performed spectroscopy using ^-polarized light and 
therefore 8 was not determined. The values for p/, and E ca ] c are taken from [28]. 


Va 

Pb 

[%] 

■Scale/ (be) 

[cm -1 ] 

■Sexp/ (be) 

[cm -1 ] 

£ 

[10 _3 cm _1 ] 

8 

[MHz] 

66 

15.60 

9388.005 

9387.9967 

8.3 

-2 

67 

28.50 

9423.589 

9423.5794 

9.6 


68 

23.21 

9454.571 

9454.5652 

5.8 


69 

7.66 

9491.049 

9491.0451 

3.9 

-2 

70 

10.82 

9525.742 

9525.7346 

7.4 


72 

39.36 

9594.454 

9594.4485 

5.5 

-22 


Table 2. Comparison of calculated (Scale) and measured (S eX p) level energies for various 
vibrational levels v/, of the Ir’ Fl„ 0 + state with / = 1, analogous to table 1. The parameter 
A is the splitting of the (A components as determined by fitting our theoretical model to the 


measured spectra (see section 6). 

Vfo 

Pb 

■E'caic/ (be) 

Eexp/ (be) 

£ 

A 


[%] 

[cm -1 ] 

[cm -1 ] 

[10 _3 cm _1 ] 

[GHz] 

73 

82.70 

9368.758 

9368.7480 

10.0 

81.84“ 

74 

69.59 

9412.519 

9412.5122 

6.8 

104.54“ 

75 

73.80 

9460.874 

9460.8718 

2.2 


76 

88.37 

9503.516 

9503.5040 

12.0 

40.01“ 

78 

58.25 

9591.479 

9591.4721 

6.9 

36.41];] 


calculation is based on a two-potential approach considering A 1 !] 3 " and b 3 Tl u (£2 = 0 + , 1,2). 
The mixing is described by the parameter pp, which represents the probability of finding the 
vibrational level in the electronic state b. Consequently, for the A state the corresponding 
parameter is given by pa — 1 — Pb■ All other admixtures like AG = 1 are negligible in our 
cases. Our spectroscopy scheme addresses only the A component of a vibrational level of 
the A — b manifold. Moreover, only states with angular momentum J — 1 and negative parity 
can be observed, because the electronic singlet component of the Feshbach molecule has the 
quantum number / = 0 and positive total parity. 
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4. Experimental observations 

4.1. Spectra of A levels 

We first discuss the data for levels with mainly A 1 ^ character. Six different vibrational states 
(v ,4 = 66 to 70 and 72) have been investigated. The obtained spectra for = 66 to 70 look 
very similar. Figure 2(a) shows the recording for va = 66 as an example. Two resonance 
dips are visible, one being the n transition (A mp — 0 ), while the other one is the o transition 
(Amp = ±1). Within the measurement uncertainty of a few MHz both resonances are located 
on top of each other and Zeeman or hyperfine splitting is not observed. We determine the 
transition frequencies from fits to the data using the function e~ KL , where the amplitude K 
is a free fitting parameter and L represents a Lorentzian. Typically, the obtained transition 
linewidths (FWHM) are on the order of 10 — 20 MHz. 

In table 1 the absolute energies of states va derived from the n resonances are 
summarized and compared to theoretical predictions. The admixing parameter pp and E ca \ c 
are taken from Ref. [28]. Since, the calculations were originally given with respect to the 
potential minimum of A 1 Eg, for the comparison to our experimental results, we added the 
electronic term energy Tf = —3993.5928(30) cm 1 / (he) of X 1 Eg [31] and the hyperfine 
shift of +8.543 GHz x h/c, where c is the speed of light. The overall agreement between 
the theoretical and experimental data is within the uncertainty of the theoretical predictions 
of 0.01 cm^ 1 /(he) (corresponding to 300MHz x h). Noticeably, the calculated values are 
systematically higher by several 10 3 cm -1 compared to our measurements. Besides a 
possible systematic uncertainty within the theoretical model, these deviations can also arise 
from the limited accuracy of our wavemeter and the uncertainty of the energy Tf. 

In contrast to the states +4 = 66 and 69, where both, the K and the a transition occur at 
the same frequency within the measurement uncertainty, +4 = 72 shows a significant splitting 
[cf. Fig. 2(b)]. This is due to the fact that the admixing of the b state is relatively large 
( p b ~ 40%, see table 1) and a b 3 H u 0 level is located energetically close-by. The level 
va = 72 significantly exhibits the characteristics of b^H u 0+, which will be discussed in the 
following sections. 

4.2. Spectra ofb levels 

Our spectroscopic data on states with mainly triplet character, i.e., pp > 50%, are shown 
in Fig. 3. For all investigated vibrational quantum numbers \p — 73 to 76 and 78 we only 
clearly observe a single resonance dip when using 7T-polarized light. Contrary to that, the 
scans related to a polarization reveal 2 or 3 resonance features of which some might have an 
unresolved substructure. In each spectrum, the a transitions are well separated from the n 
transition. We therefore choose the n resonance as a local reference to which we assign the 
frequency S — 0 in the figures. 

At first sight the spectra for \p = 73 to 76 and 78 might look somewhat irregular. For 
different vibrational levels \p the number of transitions, their splittings, and their relative 
intensities vary. In addition, the splittings for a given \p are not equidistant as mentioned 
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Figure 2. (Loss resonances for excitation of molecules from the Feshbach state to vibrational 
levels va = 66 (a) and va = 72 (b) of the A 1 Z+ potential obtained with ^-polarized light 
(black squares) and fT-polarized light (red circles). Shown is the fraction N/No of remaining 
Feshbach dimers dependent on the detuning 8 , where 8 = 0 is at the resonance frequency of 
the strong n transition. The corresponding offset energies are listed in table 1. Solid lines 
are fits of the function e~ KL to the data (see section 4.1). For a given vibrational quantum 
number va the measurements with 71 - and cr-polarized light are performed using the same 
laser intensities and pulse lengths. Colored vertical lines in part (b) indicate the frequency 
positions of the levels \F',m! p ) resulting from our model calculations (see section 6). 


earlier for high J. However, closer inspection reveals that all these spectra are characterized 
by a similar pattern. To show this, we arrange the spectra in the order v b = 74, 73, 78 and 
76 [Fig. 3(a)-(d)] corresponding to their respective splitting magnitude. In each spectrum the 
o lines are located at 8 > 0. There is always one weak resonance next to the n transition 
and one strong resonance feature at larger 8. For v b — 76, where the total splitting is very 
large, the resonance dip at 8 ~ 90MHz seems to split up into two or more lines. Due to the 
limited resolution of about 5 MHz in our experiment we cannot clearly resolve the individual 
resonance lines, but the observed fluctuations in the number of molecules are a clear indication 
of an internal structure of this resonance dip. 

In contrast, the spectrum of v b = 75 [Fig. 3(e)] is inverted compared to the spectra 
discussed before and exhibits three a resonances, all of them at 8 < 0. As it is spread 
over an even larger frequency range of about 160 MHz, the resonance dips at —160 MHz and 
— 130MHz are clearly separated from each other. The offset energies for the observed lines at 
5=0 are listed in table 2 and are compared to the theoretical predictions of Ref. [28]. Again, 
the agreement is within the theoretical uncertainty. However, we note that the measured K 
transitions contain shifts due to hyperfine and Zeeman interaction. These shifts of up to 
190 MHz (see section 6) would need to be subtracted for a proper comparison of the data 
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Figure 3. Loss spectra for excitation of molecules from the Feshbach state to vibrational 
levels Vfo = 74 (a), v/, = 73 (b), v/, = 78 (c), v* = 76 (d), and v^ = 75 (e) of the Zr 3 n„ 
potential obtained with ^-polarized light (black squares) and ct-polarized light (red circles). 
All parameter denotations, the fit function and the meanings of the vertical lines are identical 
to those of Fig. 2. The offset energies corresponding to the transitions at 8 = 0 are given in 
table 2. For v/, = 73 and 74 the intensity and pulse length of the cr-polarized light was the 
same as for K polarization. The spectra of \'f, = 78 (\'f, = 75) were measured with different 
pulse lengths T, where the ratio was T a /i n = 5/3 {T a /x n = 2/1). Concerning v^ = 76, data 
of two scans with a-polarized light are shown (magenta and red). Whereas the magenta data 
points were obtained using the same pulse area as for K polarization, it was by a factor of eight 
larger when measuring the red data points. 
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with the calculations of Ref. [1]. 

5. Simple model of the molecule 

In principle, hyperfine and Zeeman interaction within the 3 n state of diatomic molecules 
has been theoretically investigated in depth (see, e.g., the 4th-order perturbation approach 
of [33]). However, properly applying such theoretical (and often complex) approaches to 
interpret measured spectra can still be a challenge because of the large number of parameters 
for representing the different orders. Therefore, we have developed a simple model which 
neglects some fundamental properties of a molecule. Nevertheless, it should be adequate to 
explain semi-quantitatively the Zeeman and hyperfine structure observed in our spectra. 

In our model, the Rb 2 molecule is treated as a rigid rotor with fixed intemuclear 
separation. Consequently, there is no vibrational degree of freedom. However, the positions 
of the nuclei can be interchanged. This is necessary in order to construct fully antisymmetric 
wave functions for the system of the two nuclei and the two valence electrons owing to the 
particles’ fermionic character. Essentially, we consider the molecule as if it was composed of 
two unperturbed neutral atoms, of which, however, the angular momenta L and J are strongly 
coupled to the rigid rotator axis. In each of the atoms the orbital angular momentum L, of 
the local electron i is a good quantum number. Thus, molecules belonging to the atom pair 
5 Si /2 + 5 P1/2 have both a p-orbital with Li — 1 and a 5 -orbital with L\ — 0, and the total orbital 
angular momentum is L — 1 (L = L\ + Lo)■ Therefore, the two valence electrons can never 
be found in the same orbital. Coupling the electrons to the rotator axis (which corresponds to 
the intemuclear axis) forms the electronic states 2341 A i( /„ of the molecule. For simplicity, in 
the following discussion we restrict the model to those electronic states that are most relevant 
to describe our observations, i.e., states with w-symmetry and Q. = 0. These are b 2 Yl u (0 + as 
well as 0 _ ), and (2) 3 Z+ (see section 3.2). 

The molecule is described by the Hamiltonian 

H = tfoiag + //so + H r + H hf + H Z , (1) 

which, in addition to a diagonal energy matrix, contains spin-orbit coupling, nuclear rotation, 
hyperfine and Zeeman interaction. The diagonal energy matrix 

tf Diag = I E^P^+'AuO*) ( 2 ) 

A,5,0± 

sets the initial values for the energies E AS 0 ± of the electronic levels 25+1 A„ (P before the 
remaining terms of the Hamiltonian are turned on. Here, P( 2S+l A u (P) denotes the projector 
onto the respective state. In order to describe the hyperfine and Zeeman structure for a given 
vibrational level v’ (with symmetry 2S ' + 1 A!, (P), the influence of all surrounding vibronic 
levels for each symmetry is mimicked by a single, effective energy value E AS0 ±. As an 
example, let us assume that we want to describe the hyperfine and Zeeman structure of the 
vibrational level v^, = 75 of the b state. As can be seen in Fig. 1(b), v/, = 75 is surrounded by 
several v ,4 levels in its proximity, with v ,4 = 67, 68 and 69 being the closest ones. All these V 4 
levels are replaced by a single effective vibrational level with energy E^qq+ in our model. 
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The second term of Eq. (1) is the spin-orbit interaction 

Hso = Cso(Si ■ Li + S 2 ■ Li), ( 3 ) 

which couples spin 5, and orbital angular momentum L, of electron i. Here, Cso denotes the 
spin-orbit parameter being the corresponding atomic value divided by two because we have 
only 50% probability for each electron to be in the p-orbital. From the atomic fine structure in 
87 Rb (see, e.g., [34]) one obtains Cso = [E{5 2 P i / 2 ) — E(5 2 Pi/ 2 )\ /(3/? 2 ) = 2374GHz x h/h 2 . 
We use this value of Cso for the spin-orbit interaction between (2) 3 E+ and b 3 n„ 0 . These 
states are separated by about 5000 cm 1 [cf. Fig. 1(a)]. The corresponding level repulsion 
shifts the b 3 U u 0 component to lower energies by several tens of GHz x h compared to 
the situation, when spin-orbit interaction is ignored. For the spin-orbit coupling between 
A l L+ and b 3 U u 0+, we additionally take into account the overlap integral of the relevant 
vibrational wave functions, which is typically ~ 0.1 for states of the considered frequency 
range (9360 — 9600cm *). The spin-orbit interaction is responsible for the frequency splitting 
A between the 0 and 0 + components of b 3 U u and the mixing of the A and b state which is 
expressed in terms of the admixing parameter p b . It turns out that A and p b are the two 
quantities, which essentially determine the hyperfine and Zeeman structure of a vibrational 
state. By fine tuning A and p b in our model we can describe the observed spectra. For 
practical purposes, we vary neither Cso nor the value of the overlap integral (= 0.1), instead 
we use the term energies of the relevant uncoupled states in Hj) iag . Concretely, we adjust p b 
by setting the separation between A 1 E+ and b 3 Tl u , while the size of A is adjusted by shifting 
the term energy of the 0 level relative to the 0 + level. 

The third term of Eq. (1), 

H r = B w R 2 : (4) 

describes the rotation of the atom pair. According to the calculations of Drozdova et al. [28], 
the rotational constant B v is about 0.54 GHz x h/h 2 for the A — b states with dominant b 
character and vibrational quantum numbers v^ ~ 70 — 80 of 87 Rb 2 - The quantum number 
of angular momentum R appearing in the atom pair basis determines the total parity of the 
molecular state according to — (—1)^. But R is not a good quantum number for the molecular 
eigenstates since R 2 does not commute with H\y Vdg . 

Next, we consider the hyperfine interaction //hf- As mentioned in [4], the Fermi contact 
term is in general sufficient to characterize the hyperfine interaction of alkali-metal dimers. 
We use 

H HF = b F (S-T), (5) 

with S — S1+S2 and I — I\ +h- According to Refs. [4, 35, 36], the Fermi contact parameter 
b-p for an atom pair (5 + p ) is b R ~ AHF.atom/4, where Adatom = 3.417 GHz x h/h 2 denotes 
the atomic hyperfine parameter for the 5 Si/ 2 level of 87 Rb [37]. We note that this Ansatz is 
formally identical to the atomic hyperfine interaction of a ground state electron (i.e., .s' orbital) 
with its local nuclear spin The factor 1/4 normali z es the interaction because at any instant 
in time only one of the two electrons (i.e., the s electron) interacts with only one of the two 
nuclei. We fix the Fermi contact parameter to be b R = Ahf, atom/4 but note that deviations of 
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up to 20% from this approximation have been observed, e.g. for Na 2 [36]. We do not use b p as 
a free fit parameter in our model because it is strongly correlated with the unknown frequency 
splitting A (see section 6). Thus, any uncertainty in bp directly translates into an uncertainty of 
A. Furthermore, by using the ansatz of Eq. (5) we neglect the nondiagonal part of the hyperfine 
interaction with respect to S and / and thus there is no mixing of u/g symmetry. However, 
this approximation should be valid as the energy spacing between possibly coupled u/g states 
is significantly larger than the 0“/0 + spacing considered in this work. 

The last term of Eq. (1) characterizes the Zeeman interaction in a homogeneous magnetic 
field of strength B in z direction 

Hz = Mb [gL L z + gs S z + gj I Z \B , (6) 

with /1b being Bohr’s magneton. Here, we consider the Zeeman interaction due to the orbital 
angular momenta of the electrons, the electronic spins as well as the nuclear spins, where 
gL = 1, gs = 2.002319 and gj — —0.000995 for 87 Rb [37] are the corresponding y-factors. 

The matrix elements are calculated in an uncoupled atom pair basis, being a properly 
anti symmetrized product of eigenstates of all needed angular momenta and their projection 
on the space-fixed axis z, and the nuclear positions at both ends of the rotator axis. Table 3 
gives an overview of the range of quantum numbers for the Q. = 0 states of /EH,,, which are 
needed to setup the matrix. The total molecular angular momentum J (J — R + L + S, i.e., 
without nuclear spins) is a fairly good quantum number, because the hyperfine and Zeeman 
interaction is small compared to the other interactions. 

6. Model calculations and interpretation of measured data 

In the following we use the model introduced in the previous section to calculate the Zeeman 
and hyperfine structure for an A — b bound state as a function of the frequency splitting 
A = [E(0~,J = 0) — E(0 + ,J = 1 )\/h of the 0 ± components of b 3 n„:j:, the degree of mixture 
pi, between the A and b states, as well as the magnetic field B. We define A to be the splitting 
between 0^ after diagonalization of the Hamiltonian H of Eq. (1). In order to keep the 
discussion simple, we restrict ourselves to the range of quantum numbers and parameters 
directly related to our experiments. As explained in section 3.2, starting from Feshbach 
molecules we can only optically excite A — b bound levels through the A 0 + component 
with angular momentum J = 1 and negative parity. Furthermore, we want to point out that 
the eigenstates of the Hamiltonian of Eq. (1) are eigenstates of the total nuclear spin I. Thus, 
we can restrict ourselves to bound states with I — 2, being equal to the value of the singlet 
component of the Feshbach state applying the electric dipole selection rule A1 — 0. 

The diagonal Zeeman and hyperfine interactions of the states O* 1 are negligible compared 
to our measurement uncertainty. However, both the Zeeman interaction Hz and the hyperfine 
interaction //hf couple 0 + and 0 - within b 3 Tl u . In particular, (J = 1, / = 2) of 0 + couples to 

f The different J values are required for obtaining the same parity for the mixed states. In addition to J = 0, 
also the level 7 = 2 of 0~ will couple to 7 = 1 of 0 + . The rotational energy splitting between 7 = 0 and 2 of 0 
is determined by Eq. (4). 



13 


Table 3. Overview of the range of quantum numbers for states 0 + and 0 of b'’Y[ u with low 
angular momentum J. The +/— columns provide the total parity whereas I represents the 
total nuclear spin and R is the atom pair rotation. Note that the quantum numbers I and parity 
alternate with J. This behavior is also found for the even and odd values of R. Furthermore, 
the molecular rotation increases with J. 


J 

I 

0+ 

+/- 

R 

I 

0 

+/- 

R 

0 

1,3 

+ 

1 

0,2 

— 

0,2 

1 

0,2 

— 

0,2 

1,3 

+ 

1,3 

2 

1,3 

+ 

1,3 

0,2 

— 

0, 2,4 

3 

0,2 

— 

2,4 

1,3 

+ 

1,3,5 

4 

1,3 

+ 

3,5 

0,2 

T 

2, 4,6 


(J = 0, 1 — 2) and (7 = 2, 1 — 2) of 0 - (see table 3). This leads to mixing, i.e., to the creation 
of eigenstates with a net electronic magnetic moment and thus to Zeeman and hyperfine 
splittings. Interestingly, hyperfine and Zeeman interaction amplify the line splittings in a 
cooperative way because they have matrix elements for the coupling between 0 + and 0 
similar in magnitude and equal in sign. Hence, if Zeeman and hyperfine interaction are of the 
same strength their combined effect increases the line spacings not only by a factor of two but 
by a factor of four. Our experiments are indeed close to this regime for the selected magnetic 
field of about 1000 G. 

The Zeeman and hyperfine splitting crucially depends on the frequency spacing A 
between the levels 0 and 0 + of /? 3 n„. Using standard perturbation theory, the splitting 
is estimated to be proportional to (Hz + Hup) 2 /A. We recall that the spin-orbit couplings 
to (2) 3 Z+ and generate the spacing A between the levels 0 and 0 + of b 3 H u in the 
restricted Hilbert space. 

Figure 4 depicts results of our model calculations for a vibrational level of A — b with 
80% triplet (b) and 20% singlet (A) character. The frequency positions of levels | F',m' F ) are 
shown as a function of A for two magnetic fields, B — 0G and 1000 G. Although the total 
angular momentum F' is generally not a good quantum number anymore at higher magnetic 
fields, we refer to the molecular levels by the correlated value of F' at B — 0 G. In Fig. 4 only 
those six states | F',m' F ) are plotted that, according to the dipole selection rules, are accessible 
by our spectroscopy experiment. As the Feshbach molecule has F = 2, mp = 2, levels with the 
quantum numbers m' F — 2 (m' F = 1,3) can be observed via K (<y) transitions. We choose the 
state | F' — 2, m F — 2) as energy reference (i.e., 8 — 0), because it corresponds to the strong 
K resonance in Figs. 2 and 3. This allows for convenient comparison of our calculations to 
the measured spectra. Note, at B = 0G, no nip splitting can occur and therefore only three 
different curves are discernible in Fig. 4(a). 

In our calculations, A is set by adjusting the initial energy spacing between the 0 + and 
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Figure 4. Hyperfine level structure for a vibrational state v/, with pi, = 80% for two magnetic 
fields, B = 0G (a) and 1000G (b). Shown are the frequency positions 8 of the levels | F',m' F ) 
relative to the state |2,2), as a function of A. The inset in (b) gives the relative strengths K 
for the optical dipole transitions from the Feshbach state towards the levels \F'.m' F ) of v/, at 
B = 1000 G. For |2,2) we set K = 1. For convenience, we have plotted the same data in Fig. 7 
in terms of 1 /A. This makes it easier to read off the line splittings for small |A|. 
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0 components of the b state in H\y VAg [see Eq. (2)]. To a good approximation, within the 
frequency ranges considered in Fig. 4 (i.e., |A| = 16— 129 GHz) the level splittings increase 
inversely with A, just as expected from perturbation theory. This can directly be seen in Fig. 7 
of the Appendix. For B = 1000 G and a small |A| of 16 GHz, the overall spreading of the 
levels reaches more than 200MHz, whereas for a large |A| > 60GHz it is on the order of a 
few tens of MHz or less. Furthermore, the ordering of the energy levels is inverted, when 
the sign of A changes. We point out that some level spacings are smaller than the expected 
linewidths and therefore cannot be resolved in the experiment. The widths of the levels are 
mainly determined by those of the A X Y^ state (~ 12MHz) and the admixing parameter py, 
because the width of the pure /AH,, state is orders of magnitude smaller compared to the one 
of A 1 . The basic structure of the calculated levels (see Fig. 4) and the level widths let us 

expect to resolve three resonance features which agrees well with our observations shown in 
Fig. 3. 

We now want to assign the experimentally observed resonances to distinct transitions. 
Besides considering the line positions we also take into account the strength of the lines. For 
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this purpose, we calculate the dipole matrix elements A^fs,(f',»/,,) f rom the initial Feshbach 
state |FS) to the final levels \F',m' F ) of the mixed b 3 H u state. As already mentioned, we 
only have to consider the singlet component of both levels. The inset in Fig. 4(b) shows the 
relative transition strengths K = M 3 S (F , m , j/M 3 s n at B = 1000 G. We can roughly group 
the six transitions into three strong ones (by c light towards final states |3,1) and |3,3), by 
K light towards |2,2)) and three weak transitions (a: towards 11,1) and \2, 1) and k: towards 
|3,2)). Our observed spectra always exhibit only one strong line for ^-polarized light (cf. 
Fig. 3). It is therefore assigned to |2,2) and used as reference level. The calculations predict 
a second resonance for K polarization which should be about one order of magnitude weaker. 
However, we did not unambiguously observe this line since its signal is easily drowned by the 
overlaying strong a lines if the achieved light polarization is not sufficiently pure. 

In the following, the O’ transitions are discussed in detail using the results shown in 
Fig. 4(b). Although the transition towards 12,1) is weak, we should be able to clearly observe 
it, since the level 12,1) is well separated from all other levels. The three remaining transitions 
(towards 11,1), |3,1) and |3,3)), however, are quite close to each other. Especially for larger 
values of |A| (> 40GHz) these levels cannot be resolved and only a single resonance should 
be visible. Among the three transitions, the one towards |3,3) is most dominant. For low 
values of |A| (< 40GHz) the strong |3,3) line splits clearly from those corresponding to |1,1) 
and 13,1) which both barely separate. This explains why our observed spectra for c-polarized 
light in Fig. 3 exhibit at most three resonance features. At this stage we have shown that the 
experimental data can be qualitatively explained by our theoretical model and that we already 
can assign quantum numbers to the measured resonance lines. 

Now, we want to carry out a more quantitative comparison of the measured line splittings 
in Fig. 3 with the model predictions. For this, we study the dependence of the energy level 
structure on the admixing parameter pb, i.e., the percentage of the /rTl i( potential in the 
vibrational state v/,. In the simulations, we set pi ? by adjusting the term energy of the bare 
state in Eq. (2). Results for A = +59 GHz are shown in Fig. 5. To good approximation, 
within the investigated range from Pb — 65% to 95% the level frequencies depend linearly on 
Pb . This makes sense as the discussed hyperfine and Zeeman interaction only appears within 
the b 3 Tl u (£2 = 0) state. 

In order to carry out the quantitative comparison for each spectrum v/ ; , we individually 
fix the splitting of the bare A ! E+ 0 + and /rTl„ 0 + levels such that the admixing parameter pi, 
equals to its literature value [1, 28], as listed in table 2. Afterwards, we fit our model to the 
measured spectrum by adjusting a single parameter, the effective term energy of b 3 H u 0~ in 
//oiag and thus the splitting A between the 0 and 0 + components of the b state. All other 
parameters of the model are kept at the values given in section 5. For the fit, we ignore the 
states 13,2) and |1,1) since they cannot be experimentally resolved [see inset of Fig. 4(b)]. The 
resulting spectral positions of the hyperfine levels are shown in Fig. 3 as vertical lines together 
with the measured spectra. The agreement is quite satisfactory as the experimental and 
calculated line positions do not differ by more than a few MHz. Our fit results for A, i.e., the 
splitting of the states after diagonalizing the Hamiltonian H of Eq. (1), are listed in Fig. 3 
as well as in table 2. We obtain values ranging from A = —19.7 GHz to 104.5 GHz. The error 
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Figure 5. Dependence of the hyperfine and Zeeman structure on the admixing parameter /?/, 
for a fixed value of A = +59 GHz and a magnetic field of B = 1000 G. The frequency <5 is 
given relative to the |2,2) level. 


boundaries for A in table 2 are estimated by simulations shifting the resonance frequency 8 of 
the strong \F' = 3,m' F = 1) line by ±5 MHz relative to the reference 12,2). Such an approach 
is reasonable as the frequency stability in our measurements is ±(2 — 5)MHz. In addition to 
these error boundaries there exists a further uncertainty in A, since the precise value of the 
hyperfine interaction parameter hp is unknown. As already discussed in section 5, we expect 
that our adopted value for 7+ possibly deviates by up to 20% from its real value. According 
to the simple estimation for the hyperfine and Zeeman splittings (pc (Hz + //hf) 2 /A), such a 
change in 7>p leads to a similar change of |A| from the fit. In our data analysis, we have verified 
that the fits to the measured spectra remain of similar quality when we vary hp by a few tens 
of %. 

Next, we investigate the Zeeman and hyperfine structure as a function of the magnetic 
field. Figures 6(a) and (b) depict the model calculations for the vibrational levels v/ ; = 75 and 
73, respectively, using the values of pb and A given in table 2. For both plots the admixing 
parameters are similar (pb « 80%), while the respective frequency spacings A have different 
signs and magnitude. We show all the levels accessible in our spectroscopy together with the 
experimentally derived levels at B — 999.9 G. Here, the frequency reference is represented by 
the level 12,2) at this magnetic field. We present in Fig. 6(c) the full hyperfine and Zeeman 
structure of a single state v& for A = +14.0 GHz, pb = 80.00%. This graph reveals particularly 
well the transition from the linear Zeeman effect to a quadratic behavior above a few hundreds 
of gauss and the enhancement of the splitting by the cooperative effect between Zeeman and 
hyperfine interaction. 

Finally, we want to give a quantitative interpretation of the spectrum corresponding to 
the vibrational level = 72 of A 1 [see Fig. 1(b)], which has a large b state admixture 
(pb — 39.36%) due to the strong coupling of v^ = 72 to \b = 78 as these states are fairly close 
to each other. The separation is only 89.2 GHz according to tables 1 and 2. For v/, = 78, 
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B (G) 

Figure 6. Zeeman structure of the hyperfine levels \F' ,m' F ) as a function of the magnetic field 
B. The frequency is referenced to the position of level |2,2) at B = 999.9 G. (a) Simulations 
for A = -19.7 GHz, p b = 73.80%, (b) for A = +81.8 GHz, p h = 82.70% and (c) for A = 
+ 14.0GHz, pi, = 80.00%. The black square (red dot) plot symbols indicate the experimentally 
observed resonances obtained with ^-polarized (<7-polarized) light corresponding to v* = 
75 (a) and 73 (b). The error bars represent the measured transition linewidths (FWHM) 
determined from our fits to the data [cf. Fig. 3(b) and (e)]. In (c), all hyperfine energy levels 
of the vibrational state v b (J = 1, 1 = 2) are plotted, where the color code is extended for the 
additional levels compared to (a) and (b). 
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Table 4. Energy levels of the component £2 = 0 of and of b' FT,, without (w/o SO) and 
with (w SO) spin-orbit coupling. Column 3 shows the calculated energies for 0 + of the A 
and b states for 7=1, whereas column 4 contains only energy levels for 0 and 7 = 0. In 
column 5 the corresponding values for 0 + and 7 = 1 are provided, obtained with the spin- 
orbit interaction between states A and b. Column 6 reports the splittings [E(0~) — E(0 + )]/h 
calculated for the restricted system A — b and denoted by Aa-i,- while column 7 lists the results 
derived from our experiments. (For details see text.) 


state 

Va, 

0+, / = 1 

o-,/ = o 

0+,/ = 1 

Aa— b 

Aobs 

Aobs Aa~ b 



(w/o SO) 

(w/o SO) 

(w SO) 






[he x cm -1 ] 

[he x cm -1 ] 

[he x cm -1 ] 

[GHz] 

[GHz] 

[GHz] 

a 1 ^ 

65 

9352.911 


9352.078 

534.3 



b 3 n u 

73 

9369.939 

9369.901 

9368.758 

34.3 

81.8 

47-5+1°/ 

A^+ 

66 

9388.026 


9388.005 

-542.7 


b 3 n u 

74 

9414.679 

9414.641 

9412.519 

63.6 

104.5 

40-91S3 

A^+ 

67 

9423.000 


9423.589 

-268.3 



A^+ 

68 

9457.834 


9454.571 

137.9 



b 3 n u 

75 

9459.210 

9459.172 

9460.874 

-51.0 

-19.7 

- 2 + 0.6 

J 0.6 

A^+ 

69 

9492.527 


9491.049 

373.0 


b 3 n u 

76 

9503.529 

9503.492 

9503.516 

-0.7 

40.0 

40.7 til 

A^+ 

70 

9527.078 


9525.742 

655.2 



b 3 n u 

77 

9547.635 

9547.598 

9547.629 

-0.9 



AX 

71 

9561.487 


9560.041 

-373.0 



b 3 n u 

78 

9591.525 

9591.488 

9591.479 

0.3 

36.4 

36.11)4 

AX 

72 

9595.754 


9594.454 

-88.9 

-52.8 

36.11)4 

A^+ 

73 

9629.979 


9627.657 

224.9 




our model determines a spacing of A = 36.4+2 9 GHz between its 0 + and 0 components. 
Consequently, the b 3 H u 0 state is only separated by — 52.8+2 9 GHz from the A ! E+ 0 + state. 
From this, we can predict the hyperfine structure for V 4 = 72 with our model. The results are 
shown in Fig. 1(b). As can be seen, the calculated and measured resonances agree well, which 
nicely confirms the consistency of our model. 

7. Splitting between 0 + and 0~ components in a potential scheme 

In the previous section, we have determined the effective splitting A = [E(0~,J = 0) — 
i?( 0 + ,i = 1 )\/h for the state b 3 H u using our simple model without vibrational degree of 
freedom. Here, we compare the obtained results to those of coupled channel calculations with 
a potential scheme, i.e., including the full dynamics of the relative motion within the atom pair. 
As a first step we follow Ref. [1] and therefore restrict the calculations to the A — b system, 
such that the spin-orbit interaction and the molecular rotation only couple the A and b states. 
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Thus, the influence of the spin-orbit coupling of b 3 H u 0 to (2) 3 E+ 0^ is not yet considered. 
The term values for the A — b complex are listed in table 4. For the calculations we take the 
model potentials and spin-orbit functions reported in Ref. [1]. Columns 3 and 4 show the term 
values for 0 + , J = 1 and 0~, J = 0 in the absence of spin-orbit coupling, respectively. At this 
stage, the 0 + and components of the b 3 H u state are only split due to rotation. Column 5 
lists the term energies for 0 + , J — 1 if the spin-orbit coupling is included. These are the same 
values as given in tables 1 and 2. Column 6 provides the energy difference of each 0 + state in 
column 4 relative to the closest 0 state in column 3. This quantity, which we call A A _b, is the 
prediction for the 0 ± splitting within the A — b system. Noticeably, A a _ b varies strongly for 
different vibrational levels, which results from the fact, that the uncoupled vibrational ladders 
v A and v/, of A 1 Z t | and /rTl ;( are interwoven with different spacings. Thus, their repulsion 
due to the interaction of different vibrational levels is fairly irregular. Column 7, labeled A 0 b s , 
recalls the findings for A inferred from our measurements. There is disagreement between 
the measurements and the results of the restricted coupled channel model. All values A 0 b s 
are significantly larger (more positive) than A a b- However, it is striking that the differences 
between A 0 b s and A A _b, are all close to 40 GHz. The actual values are given in column 
8 , where the uncertainties are taken from A 0 b s (see, e.g., table 2). Indeed, all differences 
(except for v/ ; = 75) are equal within the given uncertainties. Could the spin-orbit coupling 
to state (2) 3 Z+ be responsible for this discrepancy of « 40 GHz? No, it cannot. Indeed, 
the (2) 3 Z+ is far up in energy [cf. Fig. 1(a)] and therefore all vibrational levels of b 3 TI u 0 
experience an almost constant shift, but it has the wrong sign to explain the observations. If 
ab initio calculations for the ( 2 ) 3 L+ potential energy curve and the corresponding spin-orbit 
interaction are applied, we obtain quantitatively that the 0 component of the b 3 U u state is 
shifted by about 50 GHz x h to lower energies. This increases the deviation of already 40 GHz 
to about 90 GHz, much beyond experimental uncertainties. 

Because the energetic order of the 0 + components are accurately determined from former 
spectroscopic work [ 1 ], our observation directs clearly to the spin-orbit energy of the 0 
component of state b 3 U u . Interestingly, the spin-orbit functions of b 3 Yl u presented in Fig. 6 (a) 
of [ 1 ] are different for Q. — 2 and = 0 by about 1.8 cm 1 as derived from the reported 
amplitude D s e ° of the spin-orbit function in Table V of Ref. [1], This difference indicates 
that b 3 U u is no longer a pure spin-orbit multiplet of a n state. If we consider that a similar 
deviation of the multiplet structure also exists for 0 ~, then a decrease of about 2.6 cm 1 
compared to the spin-orbit function for 0 + in [ 1 ] would be sufficient to model the A 0 b s as 
given in table 4. We believe that this is a plausible explanation of our observation. 

Finally, in view of the calculated splittings of 0 and 0 + we want to recall the 
measurements for the vibrational levels v A = 66 and 69 of A 'Ej (cf. Fig. 2(a) and table 1). 
The absolute values |A A _b|, even when corrected by the above determined shift of 90 GHz, are 
very large. Therefore, according to the discussion in section 6 , we do not expect significant 
hyperfine and Zeeman splittings. Indeed, no hyperfine structure was observed in our spectra 
for v A = 66 and 69. 
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8. Conclusion 

We have investigated the combined hyperfine and Zeeman structure in the spin-orbit coupled 
— b 3 Yl u complex of 87 Rb 2 dimers. We performed spectroscopy of ultracold Feshbach 
molecules at a magnetic field of 999.9 G and recorded the spectra for several excited 
vibrational levels with either dominant A l Y,+ character (v^ = 66 — 70 and 72) or dominant 
b 3 Tl u character (v/, = 73 — 76 and 78). We observe large line splittings of up to 160 MHz 
and find that the Zeeman and hyperfine structure of the 0 + state of b 3 Yl u varies strongly 
for different vibrational levels. Using a simple model, where the molecule is treated as 
a rigid rotor of two neutral atoms, we can explain the level structures as resulting from 
nondiagonal hyperfine and Zeeman interactions between the 0 + and 0 " components of b 3 H u . 
The hyperfine and Zeeman interactions act in a cooperative way, which enhances the level 
splittings. Furthermore, the level splittings depend linearly on the admixture pt> for the 
components 0 + of the complex b 3 U u and A 'E+, and scale inversely with the frequency 
spacing A between 0“ and 0 + . From fits of our model to the data, we find that A strongly varies 
in the range of —53 to 105 GHz within the interval of studied v/, and v^. Our observed values 
for A systematically deviate by about 90 GHz from predictions of close-coupled channel 
calculations of the electronic structure correlated to the atom pair asymptote 5 2 S + 5 2 P using 
empirical potentials for A ! L+ and b 3 Tl u and ab initio results for (2) 3 E+. We can eliminate 
this deviation by introducing a spin-orbit coupling function that is different by 2.6cm 1 for 
0~ and 0+. 

The fact, that we can extract the frequency spacing A from our measurements of the 
hyperfine and Zeeman structure is a significant result of this work. In ordinary spectroscopy, 
starting from a ground state molecule only the state 0+ can be addressed, owing to selection 
rules. Therefore, very little data on the A-type doubling of Q. — 0 is available, so far. 
Recently, all multiplet components of the (1) 3 IT, state were observed in the photoassociation 
of ultracold 85 Rb atoms [38]. The colliding 2 S+ 2 S atom pair contains a superposition of 
states with both symmetries 3 Z+ (0~ and 1) from which (l) 3 n ? (£2 = 0 ± , 1,2) can be reached 
by electric dipole transitions. To our knowledge, our method represents the first work to 
experimentally access A by investigating only one of the 0 ± components of b 3 H u . 

In the simple model only the hyperfine contribution from the s electron was considered. 
Thus it is unknown, how much the missing contribution by the p electron and the electric 
quadrupole interaction might influence the determination of the energy splitting between 0 
and 0+. In this respect and according to the evaluated multiplet structure here, it would 
be interesting to compare the results of the simple model to full close-coupled channel 
calculations based on accurate Born-Oppenheimer potential curves including hyperfine and 
Zeeman interaction with its non-diagonal contributions between the multiplet components. 
But for this task more experimental data is needed in order to determine the vibrational 
dependence over a large interval, which allows for deriving r-dependent functional forms 
of the spin-orbit interaction and probably also of the hyperfine interaction. The latter one was 
clearly needed in an earlier study on the triplet ground state c/ 3 Ej of Rb 2 [31]. We want to 
emphasize that the strength of the present approach is, that the modeling of all observations on 
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hyperfine and Zeeman splittings are concentrated in an effective single parameter, an average 
value of the spin-orbit function for 0 . The hyperfine parameter and the g factors are kept 
constant at their atomic values. This will be altered if a wide range of vibrational levels is 
studied. 

In general, the gained information is valuable for a detailed, fundamental understanding 
of molecular hyperfine and Zeeman level structures. Despite the already existing long 
tradition (see, e.g., the early paper by Freed [33]) there is still great interest in this field 
of research as documented by recent experimental studies, for instance with respect to 
Rb 2 [30, 39, 40]. Our results can be exploited for an optimized preparation of ultracold 
deeply-bound ground state molecules either via photoassociation [41] or a stimulated Raman 
adiabatic transfer (STIRAP) starting from Feshbach molecules. 
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Appendix 

In figure 7 we show the hyperfine level structure for the parameters of Fig. 4 but as a function 
of 1/A. Figure 7 clearly reveals the inverse dependence of the level splittings on A in the 
investigated range of | A| = 16 — 129 GHz. This representation is more convenient for reading 
off the splitting for small values of A. Furthermore, it becomes obvious that the Zeeman 
effect enhances the splittings since in (b) the frequency separation of the group of levels 
characterized by F = 1 and 3 relative to the group corresponding to F — 2 is significantly 
larger than in (a). 
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Figure 7. Hyperfine level structure for a vibrational state v/, with pi, = 80% for two magnetic 
fields, B = 0G (a) and 1000G (b). Shown are the frequency positions 5 of the levels \F',m' F ) 
relative to the state |2,2), as a function of 1 /A. 
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